Closing the gap between atomic-scale lattice deformations and continuum elasticity

Crystal lattice deformations can be described microscopically by explicitly accounting for the position of atoms or macroscopically by continuum elasticity. In this work, we report on the description of continuous elastic fields derived from an atomistic representation of crystalline structures that also include features typical of the microscopic scale. Analytic expressions for strain components are obtained from the complex amplitudes of the Fourier modes representing periodic lattice positions, which can be generally provided by atomistic modeling or experiments. The magnitude and phase of these amplitudes, together with the continuous description of strains, are able to characterize crystal rotations, lattice deformations, and dislocations. Moreover, combined with the so-called amplitude expansion of the phase-field crystal model, they provide a suitable tool for bridging microscopic to macroscopic scales. This study enables the in-depth analysis of elasticity effects for macro- and mesoscale systems taking microscopic details into account.Comment: 9 pages, 7 figures, Supporting Information availabl

Controlling the energy of defects and interfaces in the amplitude expansion of the phase-field crystal model

One of the major difficulties in employing phase field crystal (PFC) modeling and the associated amplitude (APFC) formulation is the ability to tune model parameters to match experimental quantities. In this work we address the problem of tuning the defect core and interface energies in the APFC formulation. We show that the addition of a single term to the free energy functional can be used to increase the solid-liquid interface and defect energies in a well-controlled fashion, without any major change to other features. The influence of the newly added term is explored in two-dimensional triangular and honeycomb structures as well as bcc and fcc lattices in three dimensions. In addition, a finite element method (FEM) is developed for the model that incorporates a mesh refinement scheme. The combination of the FEM and mesh refinement to simulate amplitude expansion with a new energy term provides a method of controlling microscopic features such as defect and interface energies while simultaneously delivering a coarse-grained examination of the system.Comment: 14 pages, 9 figure

A coarse-grained phase-field crystal model of plastic motion

The phase-field crystal model in its amplitude equation approximation is shown to provide an accurate description of the deformation field in defected crystalline structures, as well as of dislocation motion. We analyze in detail the elastic distortion and stress regularization at a dislocation core and show how the Burgers vector density can be directly computed from the topological singularities of the phase-field amplitudes. Distortions arising from these amplitudes are then supplemented with non-singular displacements to enforce mechanical equilibrium. This allows for the consistent separation of plastic and elastic time scales in this framework. A finite element method is introduced to solve the combined amplitude and elasticity equations, which is applied to a few prototypical configurations in two spatial dimensions for a crystal of triangular lattice symmetry: i) the stress field induced by an edge dislocation with an analysis of how the amplitude equation regularizes stresses near the dislocation core, ii) the motion of a dislocation dipole as a result of its internal interaction, and iii) the shrinkage of a rotated grain. We also compare our results with those given by other extensions of classical elasticity theory, such as strain-gradient elasticity and methods based on the smoothing of Burgers vector densities near defect cores.Comment: 16 pages, 7 figure

Commensurate-incommensurate transition and domain wall dynamics of adsorbed overlayers on a honeycomb substrate

© EPLA, 2016.We introduce an effective one-mode phase-field crystal model for studying the commensurate-incommensurate transition and domain wall dynamics of the ( √ 3×√3)R30 phase found in systems such as Xe/Pt(111), or Xe and Kr on graphite. The model allows us to study large systems where the domain walls can be separated over large macroscopic distances and at the same time incorporate the microscopic details of the domain wall structures. The resulting phase diagram shows that an intermediate stripe incommensurate phase always separates the commensurate phase from the honeycomb incommensurate phases. The energy of the domain wall crossing is investigated. We also find that near a step edge, the domain walls tend to align perpendicularly to the step edge, in agreement with recent experimental observations

Striped, honeycomb, and twisted moire patterns in surface adsorption systems with highly degenerate commensurate ground states

Atomistically thin adsorbate layers on surfaces with a lattice mismatch display complex spatial patterns and ordering due to strain-driven self-organization. In this work, a general formalism to model such ultrathin adsorption layers that properly takes into account the competition between strain and adhesion energy of the layers is presented. The model is based on the amplitude expansion of the two-dimensional phase field crystal (PFC) model, which retains atomistic length scales but allows relaxation of the layers at diffusive time scales. The specific systems considered here include cases where both the film and the adsorption potential can have either honeycomb (H) or triangular (T) symmetry. These systems include the so-called (1 × 1), (√3 × √3) R30∘, (2 × 2), (√7 × √7) R19.1∘, and other higher order states that can contain a multitude of degenerate commensurate ground states. The relevant phase diagrams for many combinations of the H and T systems are mapped out as a function of adhesion strength and misfit strain. The coarsening patterns in some of these systems is also examined. The predictions are in good agreement with existing experimental data for selected strained ultrathin adsorption layers

Honeycomb and triangular domain wall networks in heteroepitaxial systems

© 2016 Author(s).A comprehensive study is presented for the influence of misfit strain, adhesion strength, and lattice symmetry on the complex Moiré patterns that form in ultrathin films of honeycomb symmetry adsorbed on compact triangular or honeycomb substrates. The method used is based on a complex Ginzburg-Landau model of the film that incorporates elastic strain energy and dislocations. The results indicate that different symmetries of the heteroepitaxial systems lead to distinct types of domain wall networks and phase transitions among various surface Moiré patterns and superstructures. More specifically, the results show a dramatic difference between the phase diagrams that emerge when a honeycomb film is adsorbed on substrates of honeycomb versus triangular symmetry. It is also shown that in the small deformation limit, the complex Ginzburg-Landau model reduces to a two-dimensional sine-Gordon free energy form. This free energy can be solved exactly for one dimensional patterns and reveals the role of domains walls and their crossings in determining the nature of the phase diagrams

Energetics and structure of grain boundary triple junctions in graphene

© 2017 The Author(s). Grain boundary triple junctions are a key structural element in polycrystalline materials. They are involved in the formation of microstructures and can influence the mechanical and electronic properties of materials. In this work we study the structure and energetics of triple junctions in graphene using a multiscale modelling approach based on combining the phase field crystal approach with classical molecular dynamics simulations and quantum-mechanical density functional theory calculations. We focus on the atomic structure and formation energy of the triple junctions as a function of the misorientation between the adjacent grains. We find that the triple junctions in graphene consist mostly of five-fold and seven-fold carbon rings. Most importantly, in addition to positive triple junction formation energies we also find a significant number of orientations for which the formation energy is negative

Thermal conductivity decomposition in two-dimensional materials: Application to graphene

Two-dimensional materials have unusual phonon spectra due to the presence of flexural (out-of-plane) modes. Although molecular dynamics simulations have been extensively used to study heat transport in such materials, conventional formalisms treat the phonon dynamics isotropically. Here, we decompose the microscopic heat current in atomistic simulations into in-plane and out-of-plane components, corresponding to in-plane and out-of-plane phonon dynamics, respectively. This decomposition allows for direct computation of the corresponding thermal conductivity components in two-dimensional materials. We apply this decomposition to study heat transport in suspended graphene, using both equilibrium and nonequilibrium molecular dynamics simulations. We show that the flexural component is responsible for about two-thirds of the total thermal conductivity in unstrained graphene, and the acoustic flexural component is responsible for the logarithmic divergence of the conductivity when a sufficiently large tensile strain is applied

Multiscale modeling of polycrystalline graphene: A comparison of structure and defect energies of realistic samples from phase field crystal models

© 2016 American Physical Society. We extend the phase field crystal (PFC) framework to quantitative modeling of polycrystalline graphene. PFC modeling is a powerful multiscale method for finding the ground state configurations of large realistic samples that can be further used to study their mechanical, thermal, or electronic properties. By fitting to quantum-mechanical density functional theory (DFT) calculations, we show that the PFC approach is able to predict realistic formation energies and defect structures of grain boundaries. We provide an in-depth comparison of the formation energies between PFC, DFT, and molecular dynamics (MD) calculations. The DFT and MD calculations are initialized using atomic configurations extracted from PFC ground states. Finally, we use the PFC approach to explicitly construct large realistic polycrystalline samples and characterize their properties using MD relaxation to demonstrate their quality